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Xylem and phloem need to maintain steady transport rates of water and carbohydrates to 
match the exchange rates of these compounds at the leaves. A major proportion of the 
carbon and nitrogen assimilated by a tree is allocated to the construction and maintenance 
of the xylem and phloem long distance transport tissues. This proportion can be expected 
to increase with increasing tree size due to the growing transport distances between 
the assimilating tissues, i.e., leaves and fine roots, at the expense of their growth. 
We formulated whole tree level scaling relations to estimate how xylem and phloem 
volume, nitrogen content and hydraulic conductance scale with tree size, and how these 
properties are distributed along a tree height. Xylem and phloem thicknesses and nitrogen 
contents were measured within varying positions in four tree species from Southern 
Finland. Phloem volume, nitrogen amount and hydraulic conductance were found to be 
concentrated toward the branch and stem apices, in contrast to the xylem where these 
properties were more concentrated toward the tree base. All of the species under study 
demonstrated very similar trends. Total nitrogen amount allocated to xylem and phloem 
was predicted to be comparable to the nitrogen amount allocated to the leaves in small and 
medium size trees, and to increase significantly above the nitrogen content of the leaves 
in larger trees. Total volume, hydraulic conductance and nitrogen content of the xylem 
were predicted to increase faster than that of the phloem with increasing tree height in 
small trees (<~10m in height). In larger trees, xylem sapwood turnover to heartwood, if 
present, would maintain phloem conductance at the same level with xylem conductance 
with further increases in tree height. Further simulations with a previously published 
xylem-phloem transport model demonstrated that the Munch pressure flow hypothesis 
could explain phloem transport with increasing tree height even for the tallest trees. 



Keywords: carbon allocation, metabolic scaling, nitrogen allocation, phloem transport, pipe model, xylem 
transport 



INTRODUCTION 

Structural-functional tree models describe formation and growth 
of new stem and branch axes as a function of their local envi- 
ronment and topological position within a tree (e.g., Sievanen 
et al, 2000). This approach has been shown to be a powerful 
tool in reproducing realistic tree architectures (Prusinkiewicz, 
2004). Scaling rules, such as the pipe model (Shinozaki et al., 
1964), have been used to quantify secondary growth to analyse 
tree mass balance during growth (e.g., Perttunen et al., 1996). 
While these empirical rules work well to describe average tree 
growth, they actually describe fundamental functional proper- 
ties of trees related to material uptake and their long distance 
transport (e.g., West et al, 1999; McCulloh et al., 2003; Holtta 
et al, 201 1; Nikinmaa et al., 2013). Simultaneously, axial and sec- 
ondary growth patterns determine the resource allocation within 
trees. Understanding the functional implications of axial scaling 
patterns is thus an essential feature of structural-functional tree 
models. 

The long distance transport capacity of tree xylem has to be 
high enough to be able to deliver water from soil to leaves at the 
rate that it is transpired. Similarly, the long distance transport 



capacity of phloem has to match the production rate of carbon 
compounds assimilated in photosynthesis in the leaves. In addi- 
tion, water and carbon exchange rates are strongly coupled as 
both occur through the same stomatal pores in leaves, although 
water use efficiency, i.e., the ratio of plant photosynthetic produc- 
tion to plant transpiration rate, may somewhat vary across tree 
species, climatic conditions and tree size (Cernusak et al., 2007). 
It is therefore reasonable to expect that the xylem and phloem 
transport capacities, i.e., conductance per leaf area, must also be 
coupled, but they may vary within species and tree sizes. It is well- 
recognized that xylem transport capacity will ultimately limit the 
photosynthetic production rate of a tree (Tyree and Sperry, 1988; 
Jones and Sutherland, 1991). Similarly, when the transport rate of 
photosynthates through the phloem is not able to keep up with 
the rate of photosynthesis, carbohydrates will start accumulating 
in the leaves and will cause down-regulation of photosynthesis 
and/or stomatal closure (e.g., Paul and Foyer, 2001; Nikinmaa 
et al, 2013). 

The allometric relationships concerning the xylem have 
been under rigorous study for centuries [(Leonardo's notes 
MacCurdy, 2002; Huber, 1928; Zimmermann, 1983)]. Scaling 
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relationships for xylem hydraulic conductance have been devel- 
oped from theoretical basis and also tested empirically. Xylem 
sapwood cross-sectional area has been found to be conserved 
in branching, and also to be approximately linearly propor- 
tional to leaf area (Shinozaki et al., 1964; Berninger et al., 

2005) . Xylem conduit radius and conductivity has been found 
to increase from tree apex downwards, and also to increase 
with tree size (West et al, 1999; McCulloh et al, 2003) to 
compensate for the increased transport distance from soil to 
leaves. However, these changes in conductance compensate only 
partially the increased tree size (e.g., Makela' and Valentine, 

2006) . 

Allometric relations of xylem have been of great interest not 
only from the water transport point of view, but also because 
the xylem construction and maintenance consumes a major pro- 
portion of the carbon budget of a tree (Makela, 1986; Nikinmaa, 
1992). The thick lignified xylem cell walls require large amounts 
of carbon over the long transport path in trees (Hacke et al., 
2001). Xylem transport physiology and its connection with struc- 
ture are relatively well-understood. The allometric relations con- 
cerning phloem transport have received much less attention. 
Anatomical and physiological measurements are much more 
difficult to conduct on the phloem, which is a more hetero- 
geneous tissue and sensitive to external perturbations caused 
by direct measurements. Also, the transport rate of fluids, and 
therefore the required transport capacity, in the xylem typi- 
cally exceeds that of the phloem by more than an order of 
magnitude (Holtta et al., 2009). The few studies which have 
measured phloem dimensions report increased phloem alloca- 
tion and decreased conduit dimensions, i.e., cross-sectional area 
(Quilho et al., 2000), toward the tree apex. Phloem conduit 
size at tree base has also been found to increase with increases 
in tree height (Jensen et al, 2011, 2012; Mencuccini et al., 
2011). 

Although at static evaluation, the phloem comprises only a 
minor proportion of the carbon allocated to the tree, the yearly 
difference in allocation compared to xylem is not as large as 
functional phloem structures do not accumulate in the stem at 
the same extent as the xylem sapwood does. Further, phloem 
contains significant amounts of nitrogen since it is metabol- 
ically more active tissue with high protein content. Nitrogen 
is the most significant limitation of tree growth in boreal 
forests (Chapin et al., 1987; Bergh et al, 1999). Unlike total 
leaf area, which saturates to a rather stable value with increas- 
ing tree size at least at the stand level, e.g., Vanninen and 
Makela (1999), the amount of xylem and phloem tissue should 
increase with tree height. Consequently, a growing proportion 
of tree's nitrogen could be expected to be found in the trunk 
and branches of a tree with increasing tree size (Helmisaari, 
1995). 

It has been suggested that, as a growth limiting substance, 
nitrogen is allocated optimally among leaves in tree crowns 
(Mooney and Gulmon, 1979; Field, 1983). The attention has 
focused especially on optimizing photosynthetic N-use effi- 
ciency (PNUE), the ratio of photosynthesis to leaf nitrogen 
content. In the optimal solution PNUE should be maximal 
and the same for all leaves independently of their position 



in the canopy (Mooney and Gulmon, 1979; Field, 1983). Yet, 
observed nitrogen allocation to leaves typically deviates from 
theoretical expectations (Field, 1983; Hirose and Werger, 1987; 
Evans, 1993; Hollinger, 1996; Kull, 2002; Wright et al., 2004). 
One overlooked aspect that may contribute to this devia- 
tion is that while foliar photosynthetic capacity depends on 
foliar nitrogen content, balanced structural design requires that 
there has to be sufficient capacity to transport photosynthates 
from the leaves as well. One should therefore consider also 
the nitrogen required for transporting substances to and from 
the leaves (in addition to the nitrogen needed for leaf func- 
tions) to understand the optimal within tree allocation of 
nitrogen. 

We set out to determine allometric relations for the amounts 
of xylem and phloem, and their hydraulic conductances and 
nitrogen concentrations as a function of tree stem/branch diam- 
eters. We studied whether these relations were similar in four 
different species co-occurring in Southern Finland; birch (Betula 
pendula), aspen (Populus tremula), pine (Pinus sylvestris), and 
spruce (Picea abies). In general, the same biophysical constraints 
with respect to both xylem and phloem transport (in terms 
of mass flow driven by pressure gradient in a medium where 
the transport resistance is characterized by the transport dis- 
tance and the cross-section area of the conducting tissues and 
the radius of the conduits within these tissues) have to apply to 
all of these species which are living in the same environmen- 
tal conditions. We were interested whether the species specific 
differences, from e.g., differences in water use efficiency, pheno- 
logical development, nitrogen requirements, and growth rates, 
would be large enough to reflect on these allometric relations 
despite the common driving gradients. Based on the allometric 
relations derived from the measurements, we derived scaling rela- 
tions to estimate how whole tree xylem sapwood and phloem 
volume, nitrogen content and hydraulic conductance change as 
a function of tree size and how these properties are axially dis- 
tributed within a tree. We also set out to estimate how much 
nitrogen was allocated to the xylem and phloem in compar- 
ison to the leaves. Then, using a previously published xylem 
and phloem transport model, we set out to answer whether the 
Munch pressure flow hypothesis would be a consistent expla- 
nation of phloem transport of photosynthates with increasing 
tree size. 

MATERIALS AND METHODS 

DERIVATION OF SCALING RELATIONSHIPS AS A FUNCTION OF 
STEM/BRANCH DIAMETER 

We expressed stem/branch diameter (d) as a power law func- 
tion of distance from apex (x) Equation (la). Xylem and 
phloem cross-sectional areas (A x and A p ), nitrogen concen- 
trations (pjv, x and pm, p), conduit sizes (r x , r p ), and the 
fraction of xylem and phloem cross-sectional areas occu- 
pied by conduits (p c , x and p Ci p ) were expressed as power 
law functions of stem/branch diameter Equations (lb-i). The 
power function form was used since it is generally used to 
describe allometric variation in biological properties (e.g., Brown 
et al., 2004). See Table 1 for the list of symbols used in the 
manuscript. 
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Table 1 | Symbols used in the manuscript. 



UP 

'6 



A b 
Aeaf 

A P 

Ap, tot 
A x 

A x , tot 

Asw 
Ahw 

d 

Ktot, p 

Ktot, X 

*P 

fcr 

i. 
*-0 

Wp, tot 
A/x, tot 

n 

Rtot, p 



fsw, max 



Vp, tot 
^x, tot 
X 

Pc, p 

Pc, x I 
PN, phloem 
PN, bark 
PN, xylem 
Y. ot; 



Cross-sectional area of bark 
Leaf area 

Phloem (living bark) cross-sectional area 

Phloem cross-sectional area summed across all branches 

of a given diameter 

Xylem cross-sectional area 

Xylem cross-sectional area summed across all branches 

of a given diameter 

Xylem sapwood cross-sectional area 

Xylem heartwood cross-sectional area 

Branch/stem diameter 

Total phloem conductance over the whole path-length 
(over whole tree) 

Total xylem conductance over the whole path-length (over 
whole tree) 

Phloem hydraulic conductivity summed across all 
branches of a given diameter 

Xylem hydraulic conductivity summed across all branches 
of a given diameter 
Tree height 

Distance from apex at which integration of whole tree 
properties are started from 
Total amount of nitrogen in the phloem of a tree 
Total amount of nitrogen in the xylem of a tree 
Number of branches, i.e., furcations at a given height or 
stem/branch diameter 

Total phloem resistance over the whole path-length (over 
whole tree) 

Total xylem resistance over the whole path-length (over 

whole tree) 

Phloem conduit radius 

Maximum sapwood radius (in cases where heartwood 

was taken into account) 

Xylem conduit radius 

Total phloem volume in a tree 

Total xylem volume in a tree 

Distance from leaf apex 

Fraction of phloem occupied by the conducting sieve 
tubes 

Fraction of xylem occupied by the conducting sieve tubes 

Nitrogen concentration of phloem 

Nitrogen concentration of outer bark 

Nitrogen concentration of xylem 

Bases for scaling relations 

Exponents for scaling relations 



yx 



A x = ai * d\ = ai : 



A p = 0l2 * dt = 012 * 



PN, x = ot 3 * d 3 = a 3 ■■ 



PN, p = «4 * «4 = 0t 4 : 



r x = a 5 * dj = a 5 = 



yx j = aivjxf 



yx J = a 2 Y2* 2 

s\ P P 5p 
yx j 3 = a 3 Y3*3 



yx J 4 = 014Y4X4 

s\ P P 5p 
yx J = a 5 Y5% 



(la) 
(lb) 

(lc) 

(Id) 

(le) 

(If) 



jP ( A P P ; 

r p = a 6 * d\ = a6 * [yx = a 6 Y5«, 

Jp ( s\ p P i 

p c> x = ct 7 * dj = 0:7 * yyx J = 017Y7X 



Pep = a 8 



MP 



(lg) 
(lh) 

(li) 



We also took into account that xylem hydraulic conductance 
decreases radially inwards to the xylem tissue due to heartwood 
formation and/or decrease in the conduit size and connectivity 
to transpiring foliage (e.g., Melcher et al., 2003). For this, we 
assumed two extreme scenarios where (A) all of the xylem was 
conducting sapwood and (B) where conducting sapwood was 
restricted only to the outermost 2 cm of the xylem (r SWi max = 
2 cm) (e.g., Sellin, 1994). We reason that the actual amount of 
conducting sapwood, and thus xylem hydraulic conductance, 
must be in between these two extreme scenarios, and exploring 
the space between these two extremes describes how sensitive the 
scaling predictions are to changes in the radial profile of xylem 
hydraulic conductance. In addition, we assumed two extreme 
scenarios for xylem heartwood nitrogen content, where (A) heart- 
wood nitrogen concentration was the same as in sapwood and 
(B) where heartwood nitrogen content was zero. Again, the actual 
amount of nitrogen in the xylem must lie in between these two 
extreme scenarios. A literature review by Meerts (2002) done on 
71 angiosperm and 22 gymnosperm species reported that the 
nitrogen concentration of the heartwood was on average 76% of 
the nitrogen content of the sapwood, but varied a lot between 
species and studies. 

We made an assumption that the cross-sectional area of 
the xylem sapwood is conserved at branching junctions. This 
assumption is from the pipe model theory formulated originally 
by Shinozaki et al. (1964). The pipe model assumption has been 
shown to hold reasonably well for the tree species used in our 
measurements (Kaufmann and Troendle, 1981; Ilomaki et al., 
2003; Kantola and Makela, 2004; Berninger et al, 2005), and also 
to result from maximizing the carbon use efficiency of xylem 
structure (Holtta et al., 2011). Note that the original pipe model 
assumption, presented e.g., in Shinozaki et al. (1964), does not 
necessarily imply that xylem conduit radius is constant within a 
tree. We further assumed leaf area to be proportional to xylem 
cross-sectional sapwood area. This assumption is also from the 
Pipe model theory (Shinozaki et al., 1964) and is supported by 
experimental evidence (Berninger et al, 2005). Note that scal- 
ing of phloem and leaf properties are affected by the assumptions 
made about sapwood turnover to heartwood. This behavior stems 
from the pipe model assumption. The amount of heartwood 
affects the number of furcations [i.e., n in Equation (A18)] , which 
in turn affects phloem properties [see Equations (A20), (A22), 
and (A24)]. The more heartwood there is, the higher the furca- 
tion number (n) is at any given height, and the more phloem 
tissue there is. The coefficients a;, f5,-, y, and 8 in Equation (1) were 
derived from measurements or from literature estimates. We used 
the allometric relations in Equation ( 1 ) to scale whole tree xylem 
and phloem volume, conductance and nitrogen amount with tree 
height for two of the measured species (one gymnosperm and one 
angiosperm); pine and aspen. 
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TREE MEASUREMENTS 

We harvested six birch (Betula pendula), aspen (Populus tremula), 
pine (Pinus sylvestris), and spruce (Picea abies) trees from forests 
surrounding Hyytiala Forestry Field Station (61° 51' N, 24° 17' 
E, 180 m a.s.l.) and two from Ruotsinkyla research forests (60° 
22' N, 25° 00' E) so that there were two trees of each species. 
Trees varied between 6.7 and 24.9 m in height and between 4.6 
and 30.6 cm in breast height diameter and they grew in even- 
aged stands whose density varied between 8.6 and 28.3 m 2 /ha in 
basal area. The trees were harvested between May and October 
2010. In addition, two additional pine trees were harvested in 
September 2013 for some additional nitrogen content measure- 
ments. All harvested trees were healthy. The cardinal points 
were marked to the sample trees before felling. Stem diame- 
ters were measured from various relative locations within the 
stem (1.0; 2.5; 5.0; 7.5; 10.0; 15.0; 20.0; 30.0; 40.0; 50.0; 60.0; 
70.0; 80.0; 85.0, and 90.0%) in East to West and North to South 
directions. Bark thickness measurement were made from tree 
trunk at breast height and at 6 m height. We numbered all liv- 
ing branches from the base of the crown to the top of the tree 
and measured their heights from the tree base and marked their 
compass direction, length, and base diameters. Branch diam- 
eters were measured beyond noticeable basal swelling and the 
distance from tree trunk and apex were measured. The crown 
was divided into segments of height and compass directions 
and a total of 10-15 sample branches were selected from each 
tree from different heights and sides of the tree so that the 
branch size distribution was evenly represented in the sample. 
Only healthy appearing and non-damaged branches entered the 
sample. 

We measured the length of sample branches from the cut sur- 
face to branch apex and measured over and under bark diameters 
from the base of the branch. Subsequently we divided each branch 
into segments. The first segment was from the base to the first 
fork and the following segments were between subsequent fork- 
ing points. The length and one to three diameters were measured 
from each segment at 10, 50, and 90% of length along with the 
bark-less diameters. Thickness of the bark was calculated as the 
difference between under bark and over bark diameter divided 
by two. Altogether, we measured 4379 branch or stem diameters 
with and without bark between 0.9 and 276.2 mm with average 
and median diameter of 6.5 and 4.2 mm, respectively. 

For measurements of the dimensions of the living bark we cut 
85 stem disks from different heights of each tree (2.5-276.2 mm 
in diameter). The cut surfaces were sanded and scanned, and 
the thicknesses of the periderm and living bark were measured 
using a self-made image analyzer program. We used 62 sam- 
ples from stem and branch disks (3.6-276 mm in diameter) for 
nitrogen the content measurements. The bark was removed from 
xylem. The periderm of bark was then removed, and the rest of 
the bark was termed as living bark, which consisted of the pri- 
mary and secondary phloem, and vascular cambium. We dried 
the periderm and living bark samples at 60° C for 72 h and xylem 
samples for 120 h and ground all dried parts with an oscillating 
mill (Retsch MM400). The total nitrogen concentration in sam- 
ples was determined by an element analyzer (Vario MAX CN) at 
the Department of Forest Sciences at the University of Helsinki. 



Allometric equations were fitted to the measured data using 
non-linear regression between sample diameter and xylem and 
bark properties with Sigmaplot (Sigmaplot for Windows ver- 
sion 1 1.0). We used ANCOVA to compare similarities/ differences 
amongst the different species in the above regressions with 
lm function of R version 2.13.0 after the data had been In- 
transformed to correspond with the assumption of the test and 
verified the correspondence with the assumptions from residual 
plots. 

ESTIMATION OF THE NON-MEASURED SCALING RELATIONS FROM 
THE LITERATURE 

Xylem conduit radius (r x ) was estimated to scale as r x ouc 0 25 (e.g., 
West et al, 1999; Anfodillo et al., 2006). There are different esti- 
mates in the literature for the scaling of phloem conduit radius 
with tree height. We used a scaling relation of r p ax 0 25 , i.e., the 
same as for the xylem. This is an intermediate scaling between 
the scaling exponent of 0.15 reported by Mencuccini et al. (201 1), 
and 1/3 reported by Jensen et al. (2011). As distance to leaf apex 
is scaled from stem diameter, xylem and phloem conduit radius 
and conductance thus scale with stem/branch diameter, as has 
been found to be the case for xylem conduits (e.g., Zimmermann, 
1983; Olson and Rosell, 2013). We assumed that a constant frac- 
tion of the xylem and phloem cross-sectional area was conducting 
lumen volume (the rest being conduit walls, parenchyma, etc.), 
i.e., P7 and Ps (exponents for the scaling equations for p Ci x and 
Pc, p) were equal to zero. This is in agreement with a formulation 
described e.g., in Savage et al. (2010). 

SENSITIVITY ANALYSIS FOR THE SCALING RELATIONS 

To demonstrate the sensitivity of the scaling of xylem and phloem 
volume, nitrogen content, and hydraulic conductance to the val- 
ues of the measured/estimated scaling exponents we conducted 
simulations where the scaling exponents Pi, $2, P3, P4, Ps, P6, 
and 8 were varied simultaneously in random (from a linear dis- 
tribution) between 75 and 125% of the values in comparison 
to values presented in Table 2. Simultaneously, the maximum 
sapwood depth (r sw , max ) was varied between 2 and 200 cm. 

SIMULATIONS OF XYLEM AND PHLOEM TRANSPORT WITH A 
NUMERICAL MODEL 

Xylem and phloem transport and the resulting within tree pres- 
sure gradients were simulated using the scaling relations obtained 
and a previously published xylem and phloem transport model 
(Holtta et al., 2009). The model calculates xylem and phloem 
pressure and sugar concentrations and their within tree axial 
gradients in steady state. Pressure differences drive xylem and 
phloem transport, i.e., flow is proportional to pressure gradient, 
and water potential equilibrium is maintained between the xylem 
and phloem. Phloem sap viscosity was made to be sugar con- 
centration dependent. Transpiration rate, phloem loading (made 
equal to photosynthesis rate) and unloading rates and soil water 
potential were given as boundary conditions, and xylem and 
phloem hydraulic conductance and tree height were given as 
structural parameters. The position of phloem unloading could 
be varied in the transport model so that we ran simulations where 
phloem unloading was made to occur either evenly along the 
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Table 2 | Measured properties as a function of stem/branch diameter for each species from the measurements: P = B * d E , where P is the 
property under consideration, B is the base, and E is the exponent. 



Quantity Corresponding 


Species 


Intercept 


Slope 


R 2 


P 


n 


scaling coefficient 




(units) 










(intercept and slope) 














Bark cross-sectional area 


Birch 


0.45 mm 2 


1.60 


0.98 


... 


553 




Pine 


0.38 mm 2 


1.61 


0.88 


... 


2050 




Spruce 


0.52 mm 2 


1.59 


0.97 


»»« 


1371 




Aspen 


0.61 mm 


1.53 


0.97 




399 




All 


0.31 mm 2 


1.69 


0.90 


... 


4377 


Xylem cross-sectional area oti and Pi 


Birch 


0.67 mm 2 


2.01 


0.999 


... 


555 




Pine 


0.82 mm 2 


1.98 


0.999 


... 


2052 




Spruce 


0.55 mm 2 


2.05 


0.999 


»»« 


1372 




Aspen 


0.71 mm' 


2.00 


0.999 




401 




All 


0.69 mm 2 


2.00 


0.999 


»»« 


4379 


Phloem cross-sectional area 012 and P2 


Birch 


0.32 mm 2 


1.59 


0.99 


... 


10 




Pine 


0.75 mm 2 


1.27 


0.98 


,,, 


62 




Spruce 


0.0086 mm 2 


2.35 


0.93 


*** 


10 




Aspen 


0.99 mm 2 


1.3 


0.94 




37 




All 


1.97 mm 2 


1.11 


0.94 


»»« 


130 


Xylem nitrogen content 03 and ($3 


Birch 


0.11% 


0.21 


0.099 


NS 


4 




Pine 


2.02% 


-0.61 


0.86 


... 


27 




Spruce 


0.15% 


-0.021 


0.0033 


NS 


8 




Aspen 


0.68% 


—0.34 


0.80 




9 




All 


1.92% 


-0.59 


0.85 


... 


48 


Phloem nitrogen content aa, and p\t 


Birch 


7.1% 


-0.59 


0.86 


... 


12 




Pine 


1.4% 


-0.22 


0.73 


... 


28 




Spruce 


0.87% 


-0.10 


0.0387 


NS 


9 




Aspen 


22 % 


—0.79 


0.90 




13 




All 


1.48% 


-0.20 


0.38 




62 


Distance from apex y and 5 


Birch 


115mm 


1.00 


0.97 




370 




Pine 


61 mm 


1.029 


0.97 




1168 




Spruce 


85 mm 


0.98 


0.96 




1144 




Aspen 


73 mm 


1.05 


0.97 




483 




All 


87 mm 


0.98 


0.94 




3165 



Stem/branch diameter Id) is in units of millimeters (mm). Significance levels *'*p< 0.007, **P< 0.01, NS= P> 0.05. 



phloem transport pathway or exclusively in the roots. We did 
three simulations with the model. (1) We used the model to 
calculate the axial xylem and phloem pressure and sugar con- 
centration gradients taking into account the axial distribution 
of xylem and phloem tissue and their specific conductivity. We 
used a 10 m pine with a maximum sapwood depth of 2 cm as 
an example, and took the axial distribution of xylem and phloem 
conductivity from the equations shown in the Appendix A3 and 
demonstrated in Figure 7. (2) We modeled how phloem tissue 
should be distributed axially in order to minimize turgor pres- 
sure difference between the leaves and roots, i.e., source and 
sink, for a fixed phloem volume at the whole tree level. In other 
words, the total amount of phloem tissue was preserved, but 
was distributed unevenly as a function of axial position. (3) We 



simulated how phloem turgor pressure difference between the 
leaves and roots would change as a function of increasing tree 
height, and whether phloem transport would be able to func- 
tion according to the Munch pressure flow hypothesis when trees 
become taller. For this, we used the scaling relations for whole tree 
xylem and phloem hydraulic conductance as a function of tree 
height derived in Appendixes Al and A2 (and demonstrated in 
Figure 4 and Table 4) for the case of pine. In this simulation, leaf 
gas exchange rates (i.e., transpiration and photosynthesis rates) 
were determined for each tree height so that leaf xylem water 
potential always remained at a constant value, i.e., we assumed 
isohydric behavior. Leaf water potential was held at — 2.0MPa, 
which is a typical value for Scots Pine in Hyytiala Forestry Field 
Station in summer conditions (Martinez- Vilalta et al., 2009). 
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Photosynthesis rate was made to be proportional to transpiration 
rate. Water use efficiency was set to 250, a typical value for Scots 
Pine in Hyytiala Forestry Field Station (e.g., Hari and Makela, 
2003). The absolute value for phloem conductance was chosen so 
that the turgor pressure difference between the leaves and roots 
obtained a reasonable value, ~0.7 MPa for the case of unloading 
in the soil for a 10 m tree. Note that this resulted in different ini- 
tial values for phloem conductance between the cases where no 
heartwood was assumed and the case where the maximum sap- 
wood depth was set at 2 cm. We also varied the initial value of 
phloem conductive and distribution of phloem unloading to see 
their effects on the results. 

RESULTS 

SCALING OF WHOLE TREE PROPERTIES FROM ALLOMETRIC 
EQUATIONS 

The equations obtained for the scaling of xylem and phloem 
properties as a function of tree height (L) (starting from a distance 
Lo from leaf apex) are as follows 
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where n is the number of branches, i.e., furcations, at a distance x 
from the leaf apex, V Xj tot is the total xylem volume in a tree, V x< to t 
is the total xylem volume in a tree, JV Xj tot is the total amount 
of nitrogen in the xylem in a tree, Nn tot is the total amount of 
nitrogen in the phloem in a tree, Nleaf, tot is the total amount of 



nitrogen in the leaves of a tree, K Xt tot is the whole tree hydraulic 
conductance of the xylem, K x ^ tot is the whole tree hydraulic 
conductance of the phloem. The derivation of the equations is 
presented in the Appendix Al. These equations apply only to the 
case without heartwood. The numerical equations for the whole 
tree scaling relations including sapwood to heartwood turnover 
are shown in Appendix A2. Lo was set to 0.1 m, except in the case 
of aspen phloem nitrogen content, in which Lo was set to 1.0 m. 
The values of 0.1 m and 1 m for Lo correspond to branch diam- 
eter of ~2 and 20 mm, respectively (see Figure 2C). Xylem and 
phloem properties were given constant values at branches than 
smaller than this. The values for Lo were chosen large enough 
so that we had measurements from branches of corresponding 
diameter. 

TREE MEASUREMENTS 
Whole bark vs. diameter 

The cross-sectional area of the whole bark {Ay), i.e., inner plus 
outer bark, increased with stem/branch diameter following a rela- 
tionship Ai, = 0.31d L68 when all the data was pooled together 
(Figure 1A). The scaling exponent ranged from 1.5344 in aspen 
to 1.59-1.61 in spruce, birch and pine with all relations being 
highly significant (see Table 2). The scaling exponents were rather 
close to each other across the species. When testing the difference, 
the logarithmic transformation changed the exponents somewhat 
(1.45, 1.50, 1.52, and 1.60 in pine, birch, aspen, and spruce, 



A 1e+5 




1 10 100 

Stem/branch diameter (mm) 



1000 



FIGURE 1 | Measured bark cross-sectional (A) and the ratio of bark to 
xylem cross-sectional area (B) as a function of stem/branch diameter. 
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FIGURE 2 | Measured inner bark (i.e., phloem) cross-sectional area (A), 
xylem cross-sectional area (B), and distance from leaf apex (C) as a 
function of stem/branch diameter. 



respectively). Out of these values, only the spruce exponent was 
very significantly different (p < 0.001) from the birch exponent. 
The ratio of whole bark to xylem diameter showed a consis- 
tent trend of decreasing with increases in stem/branch diameter, 
following a relationship Ai,/A x = 3.32d~ 094 (Figure IB). The 
relation of whole bark diameter against stem/branch diameter 
was not used in the scaling predictions since it does no. Instead, 
the bark was divided into outer and inner bark, and the latter 
represents the functional phloem tissue. 

Inner bark and xylem cross-sectional area vs. diameter and 
stem taper 

For inter-species comparison of inner bark thickness there was 
sufficient data for aspen and pine. Their exponents were not sig- 
nificantly different from each other in the ln-transferred data. 
Pooling all species data together for the non-linear regression 
yielded A p = 1. 2d 1 - 286 for the relation between bark and 



branch/stem diameter (Figure 2A, Table 2). The cross-sectional 
area of the xylem (A x ) increased with stem/branch diameter 
in slightly different manner between the species (see Table 2). 
Pooling the data together yields A x = 0.68d 2 009 for the rela- 
tion between xylem tissue and branch/stem diameter (Table 2, 
Figure 2B). When all data was pooled together, the cross-sectional 
areas of xylem and phloem were equal at branch diameters 
of ~2.5 mm. Phloem thickness resulted in saturating to a very 
constant value (phloem thickness = 0.88 *d° M2 , all species com- 
bined together, not shown). When data from the main stems plus 
the first order branches was pooled together, the distance from 
leaf apex was found to scale with stem/branch diameter following 
ad 098 (Figure 2C). 

Nitrogen content 

Nitrogen content increased clearly with decreasing stem diam- 
eter in both the living bark and the whole bark, but remained 
fairly constant for the xylem (Figure 3). When all species were 
pooled together relations of pjv, phloem = 0.0148 d~ 0 - 20 for the 
inner bark, and pjv, xylem = 0.0192 d~ 059 for xylem (Table 2). 
As there were only four birch values that differed clearly from 
the rest we also calculated the relationship without birch yield- 
ing pjv. xyiem = 0.0021 d~ 0 - 64 . While all species seemed to follow 
similar pattern for the living bark, there seemed to be a level dif- 
ference for the whole bark so that there was the most nitrogen 
in the aspen bark and least in pine bark for the same diame- 
ter. The phloem nitrogen concentrations in the smallest branch 
diameters were found to be very close to the nitrogen concentra- 
tion found in the foliage of the corresponding species in pine and 
aspen (compare Figure 2A to i ea f in Table 3). 

Whole tree scaling relations 

Whole tree scaling relation predictions were made for two exam- 
ple species: pine and aspen. The allometric relations used in 
the scaling of whole tree xylem and phloem volume, nitrogen 
content and hydraulic conductance are presented in Table 3. 
Some of the parameter values (ai_4,Pi_5, y, and 8) used in 
the scaling were taken from the measurements, while the rest 
(a 5 - a 8> p 6 - p 9 , C; s , psla, and p Ni i eaf ) were estimated from 
literature. Figure 4 shows the scaling relations for phloem and 
leaf properties in relation to the xylem properties, and Figure 5 
shows the absolute values for xylem, phloem and leaf properties. 
Whole tree phloem volume and hydraulic conductance decreased 
more sharply in comparison to whole tree xylem sapwood vol- 
ume and hydraulic conductance with increases in tree height in 
small tree heights (<10m in height), and also in larger trees 
in the absence of xylem heartwood formation (Figures 4A,B, 
5A-D,G,H)- However, at trees larger than ~10 m in height, xylem 
heartwood formation according to the maximum 2 cm sapwood 
depth scenario was found to maintain phloem conductance and 
volume at approximately a constant proportion of the xylem sap- 
wood with increases in tree height. Aspen had a larger amount 
of phloem and higher phloem to xylem ratio in relation to pine. 
Leaves were the largest sink of nitrogen in small trees, but xylem 
and phloem exceeded the leaves as a nitrogen sink with increases 
in tree height (Figures 4C,D, 5E,F). The total nitrogen content of 
the phloem was smaller than that of the xylem in pine and large 
aspen trees. The total nitrogen content of the phloem exceeded the 
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FIGURE 3 | Nitrogen concentration of living bark (A), xylem (B), and 
whole bark (C) as a function of stem/branch diameter. The data 
measurement points in living bark and xylem were from 8 trees (4 species). 



xylem nitrogen content in small aspen trees (Figures 4C,D, 5E,F). 
Assumptions on heartwood proportions and nitrogen content of 
the heartwood caused the relative nitrogen contents between the 
tissues to vary strongly. When heartwood nitrogen concentration 
was low ("hw N 0%" in Figure 4C), phloem remained an impor- 
tant nitrogen sink relative to xylem even in large trees. However, 
when there was no heartwood, or the nitrogen content of heart- 
wood was assumed to be same as that of the sapwood, then the 
role of the phloem as a nitrogen sink decreased in relation to 
xylem with increases in tree size. Table 4 present the absolute 
values for scaling of tree xylem and phloem volume, nitrogen 
content, conductance, and leaf area-specific conductance as a 
function of tree size. Note that scaling is not strictly allometric 
[see Equation (2) and Appendix A2] , although very close to it, for 
each case. 



Sensitivity analysis for the scaling relations 

Figure 6 shows the minimum and maximum xylem and phloem 
volume, nitrogen content and conductance in relation to a 
10 m tree obtained from the sensitivity analysis done with 1000 
parameter combinations. The general trends within remained 
unchanged, although the xylem, phloem and leaf properties over- 
lapped with each other. Xylem and leaf properties seemed to be 
more sensitive to parameter combination than those of phloem. 

Axial distribution within a tree 

Within a 10 m tree (taken as an example here) phloem cross- 
section (and volume) was distributed very much toward the apex, 
whereas xylem sapwood cross-section was evenly distributed axi- 
ally, following from our pipe model assumption (Figure 7A). 
Xylem and phloem nitrogen content were more concentrated 
toward the apex (Figure 7B), but this relation was much stronger 
for the phloem, especially for aspen. Xylem conductivity was 
more concentrated (k x a x 03 ) toward the base, whereas phloem 
conductivity was concentrated toward the apex when no heart- 
wood formation was assumed (Figure 7C). Assuming maximum 
sapwood depth to be 2 cm caused phloem conductance to be 
distributed more evenly within the transport axis. The axial dis- 
tribution of xylem and phloem properties was very similar in pine 
and aspen for cross-sectional area and conductance, but differed 
greatly for nitrogen content. 

Simulations for the pressure and sugar gradients within the xylem 
and phloem, simulation 1 

The xylem pressure (water potential) drop was predicted to occur 
more steeply close to the apex, while phloem pressure drop was 
predicted to occur more at the tree base in (Figures 8A,B), par- 
ticularly when phloem unloading occurred in the soil. Phloem 
pressure gradients were sensitive to heartwood assumptions. In 
the absence of heartwood formation, phloem hydraulic conduc- 
tivity was more concentrated toward the apex (see Figure 7C), 
which resulted in the phloem turgor pressure drop to concentrate 
more toward the base of the tree (Figure 8B). Phloem osmotic 
concentration gradient, which results from the interplay between 
both xylem and phloem transport properties, was predicted to 
be more evenly distributed over the transport axes. The normal- 
ized pressure and concentration gradients shown in the figure 
were not very sensitive to parameterization of the model, but the 
absolute values naturally were (not shown). 

Simulations for "optimal" allocation of axial phloem transport 
capacity, simulation 2 

When all of the sugar unloading occurred in the roots, the lowest 
turgor pressure difference between source and sink was obtained 
when phloem cross-sectional area scaled by A p a x~ 0A7 (x is dis- 
tance from apex). When phloem unloading was made to occur 
evenly within the stem, a scaling relation of A p a x~°- 33 yielded the 
lowest turgor pressure difference. Importantly, the optimal axial 
allocation of phloem tissue predicted by the model was never as 
large as in the scaling results from the measurements, i.e., A p a 
x -o.69f or p me an j ^ a x -o.67f or aS p en f or a 10 m tree with- 
out heartwood and A p a x~ 0A5 for both pine and aspen when a 
maximum sapwood depth of 2 cm was assumed. 
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Table 3 | The scaling coefficients of xylem and phloem properties as a function of stem/branch diameter used in the scaling Equation (1). 



Scaling coefficient 


Symbol 


Value used 
for pine 


Value used for 
aspen 


Measured/estimated 
from literature* 


Base for stem/branch diameter as a function of 
distance from leaf apex (mm) 


y 


15.8 mm 


15.6mm 


Measured 


Base for xylem cross-sectional area as a 
function of stem/branch diameter 


<*i 


0.82 


0.71 


Measured 


Base for phloem cross-sectional area as a 
function of stem/branch diameter 


(12 


0.75 


0.99 


Measured 


Base for xylem nitrogen concentration as a 
function of stem/branch diameter 


<*3 


2.02 


0.68 


Measured 


Base for phloem nitrogen concentration as a 
function of stem/branch diameter 


0(4 


1.4 


22.1 


Measured 


Base for xylem conduit radius as a function of 
stem/branch diameter 


<*5 


1 unit less 


1 unit less 


Irrelevant as only relative values are shown 


Base for phloem conduit radius as a function 
of stem/branch diameter 


06 


1 unit less 


1 unit less 


Irrelevant as only relative values are shown 


Raco fnr v\/lom onnrli nt Honcitv/ /ratin nf f^nnHi nt 
Uaot! I Ul Ay Id I 1 IjUI IUUIL UcllolLy \lallU Ul L.UI IU U 1 L 

cross-sectional area to tissue cross-section 
area) as a function of stem/branch diameter 


a 7 


1 1 1 n it locc 

I U I 1 1 L I Goo 


1 i in it locc 
I U I 1 1 L I Coo 


Irrelevant ac nnl\/ relatives valiioc aro chn\n/n 
n icicvai il do uiny iciaLivc vaiuco ait; oi luvv 1 1 


Base for phloem conduit density (ratio of 
conduit cross-sectional area to tissue 
cross-section area) as a function of 
stem/branch diameter 


a s 


1 unit less 


1 unit less 


Irrelevant as only relative values are shown 


Exponent for stem/branch diameter as a 
function of distance from leaf apex 


5 


0.97 


0.97 


Measured 


Exponent for xylem cross-sectional area as a 
function of stem/branch diameter 


Pi 


1.98 


2.00 


Measured 


Exponent for phloem cross-sectional area as a 
function of stem/branch diameter 


P 2 


1.27 


1.31 


Measured 


Exponent for xylem nitrogen concentration as 
a function of stem/branch diameter 


p 3 


-0.61 


-0.34 


Measured 


Exponent for phloem nitrogen concentration 
as a function of stem/branch diameter 




-0.22 


-0.80 


Measured 


Exponent for xylem conduit radius as a 
function of stem/branch diameter 


P 5 


0.25 


0.25 


West et al., 1999 


Exponent for phloem conduit radius as a 
function of stem/branch diameter 


Pe 


0.25 


0.25 


West et al., 1999; Jensen et al., 2011; 
Mencuccini et al., 2011 


Exponent for xylem conduit density (ratio of 

ffiniHi nt rrncc-cortinna 1 area tn ticci id 
UUIIUUIL LI Uoo OCLUUI lal alca LU LlooUC 

cross-section area)as a function of 
stem/branch diameter 


P 7 


0 


0 


Savage et al., 2010 


Exponent for phloem conduit density (ratio of 
conduit cross-sectional area to tissue 
cross-section area) as a function of 
stem/branch diameter 


Ps 


0 


0 


Assumed to be the same as for the xylem 


Leaf to sapwood area ratio C; s 


C/s 


3400 m 2 /m 2 


2700 m 2 /m 2 


Hoffmann and Usoltsev, 2002; 
Martinez-Vilalta et al., 2009 


Specific leaf area 


PSLA 


300 g/m 2 


80 g/m 2 


Niinemets, 1997; Martinez-Vilalta etal., 
2009 


Leaf nitrogen content 


PN, leaf 


1.2% 


2.4% 


Niinemets, 1997; Martinez-Vilalta etal., 
2009 



Stem/branch diameter Id) is in units of millimeters (mm). 
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FIGURE 4 | The predictions for the whole tree phloem volume in 
relation to xylem sapwood volume (A), phloem hydraulic conductance 
in relation xylem hydraulic conductance (B), total phloem and leaf 
nitrogen content in relation to xylem hydraulic content for the 
scenarios in which the heartwood has the same nitrogen content as 
the sapwood and for the case of no heartwood (C), and phloem and 
leaf nitrogen content in relation to xylem hydraulic content for the 
case where the heartwood has the same nitrogen content as the 
sapwood (D). "sw = 2 cm" refers to a scenario where maximum sapwood 
depth was assumed to be 2 cm and the heartwood did not contain any 
nitrogen, "no hw" to a scenario of no heartwood formation, i.e., all of the 
xylem was assumed to be conducting, and "hw full N" a scenario where 
maximum sapwood depth was assumed to be 2 cm and the heartwood 
nitrogen content was the same as that of the sapwood. In (B) the same 
area-specific conductivity was assumed for xylem and phloem. 



Simulated whole tree phloem turgor pressure differences, 
simulation 3 

Finally, in simulation 3, we analyzed how the whole tree level tur- 
gor pressure difference varies as a function of tree height using 
the predicted structural scaling of whole tree xylem and phloem 
hydraulic conductance. Phloem turgor pressure was predicted 
generally to increase slightly with increases in tree height when no 
heartwood formation was assumed, and to decrease slightly when 
maximum sapwood depth was limited to 2 cm (Figure 9). As the 
actual amount of sapwood can be predicted to lie in between these 
extreme scenarios, the turgor pressure differences between the 
leaves and roots could thus be expected to remain rather stable 
with increases in tree height. Phloem became unable to trans- 
port all of the assimilated sugars in trees larger than 15 m only in 
the case of low initial phloem conductivity and the assumption of 
no heartwood formation. In this case phloem sap viscosity expe- 
rienced a sharp build up preventing an increase in the phloem 
transport despite an increase in the turgor pressure gradient. The 
increase in turgor pressure difference with increasing tree size 
was more pronounced when sugar unloading occurred exclusively 
in the roots (Figure 9A) in comparison to phloem unloading 
occurring evenly along the stem (Figure 9B). In many of the 
cases presented, phloem turgor pressure difference increased with 
increasing tree height for small trees, but then started to decline 
again. This was due to gravity which started become important 
for taller tree. Gravity aids phloem transport while decreasing the 
capacity of the xylem to transport water to the leaves. According 
to the isohydric scenario presented here, the decrease in xylem 
transport led to lower leaf exchange rates and thus also for a 
smaller transport need for the phloem. Increase in the initial value 
for phloem conductivity decreased the turgor pressure gradient 
for all tree sizes, as would be expected. Importantly, the turgor 
pressure difference between the leaves and roots required to drive 
the phloem transport of the assimilated sugars was predicted not 
to increase linearly with increases in tree height. 

Pine was used as an example species in all of the simulations 
done with the xylem and phloem transport model, but the corre- 
sponding simulations for at least aspen would yield similar results 
as the scaling relations for the xylem and phloem volumes and 
hydraulic conductances are quite similar amongst the species (see 
Figure 4 and Table 4). 

DISCUSSION 

The equations constructed in this study make it possible to esti- 
mate whole tree level xylem and phloem properties (volume, 
hydraulic conductivity, nitrogen content). The equations require 
relations for xylem and phloem thicknesses/cross-sectional areas, 
conduit sizes and densities, and the nitrogen content across 
stem/branch samples of different diameters, and on the shape of 
the stem taper. Predictions can be made on how whole tree level 
properties scale with tree size assuming that the measured rela- 
tionships do not change with tree height. This was supported by 
the data presented here on trees that varied in size measured for 
four different species. Phloem, bark, and xylem cross-sectional 
areas, and nitrogen content could be predicted with minimal 
deviation from stem/branch diameter alone using allometric, i.e., 
power law, relationships. The approach presented here can also 
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FIGURE 5 | The predictions for the absolute values for whole tree 
volume of xylem and phloem (A,B), hydraulic conductance of xylem and 
phloem (C,D), nitrogen content of xylem, phloem and leaves (E,F), and 
hydraulic conductance of xylem and phloem per leaf area (G,H) as a 
function of tree height, "sw = 2 cm" refers to a scenario where maximum 



sapwood depth was assumed to be 2 cm and the heartwood did not contain 
any nitrogen, "no hw" to a scenario of no heartwood formation, i.e., all of the 
xylem was assumed to be conducting, and "hw full N" a scenario where 
maximum sapwood depth was assumed to be 2cm and the heartwood 
nitrogen content was the same as that of the sapwood. 



be connected to functional-structural tree models that often pro- 
vide detailed description of tree axes and their dimensions (e.g., 
Sievanen et al., 2000). 

Phloem volume and nitrogen content were predicted to be 
concentrated heavily toward the tree apex, in contrast to the 
xylem, whose properties were more evenly distributed within a 
tree (Figure 7). Partially the latter was due to the pipe model 
assumption for the xylem. However, the pipe model assumption 
has been shown to work quite well for all the species analyzed in 
the measurements (Kaufmann and Troendle, 1981; Ilomaki et al., 
2003; Kantola and Makela, 2004; Berninger et al, 2005). Also 
phloem transport capacity (hydraulic conductance) was concen- 
trated more toward the apex, especially if heartwood formation 



was limited. In contrast, xylem conductance was concentrated 
toward the base. In both cases the translocation capacity is thus 
largest closest to the source of the principal transported sub- 
stance. Our measurements on the axial profile of the phloem 
thickness and cross-sectional area confirm the findings of Quilho 
et al. (2000) who found that the thickness of conducting phloem 
cells is approximately constant along the stem height, while the 
cross-sectional area of water conducting xylem increases clearly 
from stem top toward base. The highly uneven axial distribution 
of nitrogen in both the xylem and phloem (Figure 3) signifies that 
nitrogen sampling must be done from varying branch sizes and/or 
tree heights in order to get an appropriate estimate for whole 
tree level xylem and phloem nitrogen content. For example, if 
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Table 4 | The results for scaling of tree properties as a function of tree height (L). 



Property 


Pine xylem 


Pine phloem 


Aspen xylem 


Aspen phloem 


Volume (no hw) 


T*m-5* / 2 93 

2. *1(J L_ ' 


c r-#-ir,— 5* i 2 27 


+ r\—b* i 2 94 
Z IU /_ 


n r\*-ir\ — 5 * / 2 31 

y.u IU L 


Volume (sw = 2 cm) 


n*-in-5* / 1 99t 


a a * ir\ — B* i 2 11 
4.4 IU L 


o*m-5* / 1 99t 

y* IU L 


e? i*m — 5 * / 2 14 
O.Z 10 u */_^''^ 


Nitrogen content (no hw) 


n "7Q* / 2 3467 

u. /y L 


n oc* / 2 11 
U.oO L 


n o /i # / 2 61 


1 n/i oo# / 1 9954 
1 .U4oo L -^ JiJ ^ 


Nitrogen content (sw = 2 cm, hw full N) 


n "7Q* / 2 3467 

u. /y 


n on* / 1 -88 
U.oO L 


n o /i # / 2 61 


1 no* / 1.82 
I .Uo L 


Nitrogen content (sw = 2 cm, hw N 0%) 


n rr o * / 1 F>6 
Z.OZ L 


n on* / 1 88 
U.oO L 


1 c/i * ; 1 70 
1 ,54 L 


1 no* ; 1 82 
1 .(Jo L 


Conductance (sw = 2 cm) 


a /_ 054 


a /_0.56 


a i_0.55 


a |_0.53 


Conductance (no hw) 


a/. 1 - 45 
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Only the exponents were estimated for the conductances, "(sw = 2cm, hw N 0%)" refers to a scenario where maximum sapwood depth was assumed 2cm 
and heartwood was assumed not to contain any nitrogen, "sw = 2cm, hw full N" a scenario where maximum sapwood depth was assumed to be 2cm and the 
heartwood nitrogen content was the same as that of the sapwood, and "no hw" refers to a scenario of no heartwood formation, i.e., all of the xylem was assumed 
to be conducting. 
1 Sapwood volume 



nitrogen content sampling was done exclusively from larger stem 
and branch parts, then the total amount of nitrogen allocated to 
the vascular tissues would be grossly underestimated. This result 
has also direct implications to forest management where bioen- 
ergy harvesting is becoming more popular with the need for 
boosting the use of renewable energy sources. Our results imply 
that removal of distal parts of the crown from the growing site will 
deplete the ecosystem nitrogen pool as efficiently as the removal 
of leaves. 

The relations between xylem sapwood and phloem volumes 
and conductances at the whole tree level were found to be 
sensitive to the assumption made about sapwood turnover to 
heartwood. When no heartwood formation was assumed, whole 
phloem conductance could not keep up with xylem conductance 
with increase in tree height. However, when a maximum sapwood 
radius of 2 cm was assumed, whole tree xylem and phloem con- 
ductances were predicted to change at approximately the same 
rates with tree growth, and xylem sapwood to phloem ratio was 
predicted to saturate approximately to a value of 10 (Figure 4A). 
This is approximately the same ratio as Holtta et al. (2009) pre- 
dicted based on the ratio of typical water to CO2 exchange rate 
the leaves and ranges of xylem and phloem pressure gradients 
typically observed. 

It seems clear that the xylem and phloem become increas- 
ingly larger sinks of nitrogen in relation to foliage with increases 
in tree height, and that the nitrogen requirements of the vas- 
cular tissues could be a major limiting factor to tree growth in 
the Boreal region. In trees larger than ~10m (Figures 4C,D; the 
exact height was depend on the species and the assumptions 
made about heartwood), xylem had the highest nitrogen con- 
tent. Also some previous studies have reported large amounts 
of nitrogen in the stems of large trees (e.g., Helmisaari, 1995). 
Aspen had a larger proportion of nitrogen in the phloem in 
comparison to xylem and leaf than pine. Also the proportion of 
nitrogen in the xylem in comparison to the leaves was smaller in 
aspen in relation to pine. The case of nitrogen allocation between 



phloem and foliage is particularly interesting as there is a clear 
tradeoff between the nitrogen used to assimilation or assimilate 
transport. Already Mooney and Gulmon (1979) and Field (1983) 
suggested on theoretical grounds that optimal nitrogen allocation 
within tree crowns should yield constant photosynthetic nitrogen 
use efficiency. However, such distribution has rarely been found, 
probably owing to various other factors that influence photosyn- 
thetic production rate of foliage in the crown apart from nitrogen 
(e.g., Posada et al., 2012). In reality, the proportion of nitrogen 
allocated to the leaves could decline even more strongly with 
height than our analysis suggests as the leaf to sapwood ratio 
typically decreases with increases in tree height (e.g., McDowell 
et al, 2002a,b), which was not taken into account in our analysis. 
However, not all of the nitrogen found in the xylem and phloem 
is necessarily bound to the tissue structure, but it could also be in 
temporary storage there (Wetzel et al., 1989). Our study was con- 
ducted in the boreal environment where soil water availability is 
hardly ever restricting tree function and growth, while nitrogen is 
the main resource limitation for tree growth. It would be interest- 
ing to see if the allometric relations observed here diverge if trees 
from different environments would be added to the comparison. 

The phloem transport capacity was predicted to decline more 
strongly than the xylem transport capacity when a tree grows in 
height, although the scenario of rapid xylem sapwood to heart- 
wood turnover led to the ratio of phloem to xylem to phloem to 
stabilize at tree heights larger than 10 m. Theoretically, xylem and 
phloem transport conductances should scale almost equally with 
growth in height, if the ratios between water and carbon exchange 
and the driving forces for xylem and phloem transport are to 
be maintained. How is it then possible that leaf specific phloem 
transport capacity will decrease more in proportion to xylem 
transport capacity? We can hypothesize several explanations for 
this; (A) Gravity will increase the flow rate in the phloem and 
decrease the flow rate in the xylem for a given pressure gradient 
with increasing tree height. (B) A large proportion of the photo- 
synthates might be consumed close to the apex in tall trees, so that 
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FIGURE 6 | The minimum and maximum xylem sapwood and phloem 
volume (A), nitrogen content (B) and conductance (C) in relation to a 
10 m tree obtained from the sensitivity analysis done with 1000 
parameter combinations. Also total leaf nitrogen content is shown in (B). 



phloem conductance can be allowed to decline at lower heights 
in the tree. This is in contrast to the xylem where practically all 
of the water is transported all the way from soil to the foliage. 
(C) Trees compensate for the decreased leaf area specific phloem 
conductance by increasing the turgor and osmotic pressure gra- 
dient in the phloem as a tree grows in height. Our simulations 
with the coupled xylem and phloem transport model revealed 
that the factors given above could explain the functioning of 
phloem transport by the Munch pressure flow hypothesis with 
increasing tree height even when phloem conductance decreases 
slightly in relation to xylem conductance with increasing tree 
height. The simulations also revealed that in most of the scenar- 
ios explored the turgor pressure difference between the leaves and 
roots remained rather constant with increases in tree height. This 
result is line with a recent review (Turgeon, 2010) which stated 
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FIGURE 7 | The axial distribution of xylem sapwood and phloem 
cross-sectional area (A), nitrogen content (B) and hydraulic 
conductivity in a 10-m tree (C). Values are expressed in relation to tree 
base in each case. "sw=2cm" refers to a scenario where maximum 
sapwood depth was assumed to be 2cm, and "no hw" refers to a scenario 
of no heartwood formation, i.e., all of the xylem was assumed to be 
conducting. 



that there are indications that the turgor pressure differences 
between the sources and sink would not increase with tree size. 

The xylem and phloem transport model also predicted that 
concentrating phloem volume more toward the leaf apex yielded 
lower turgor pressure difference between leaves and roots, espe- 
cially if part of the sugars transported in the phloem are utilized 
along the stem. However, the actual increase in phloem volume 
toward the apex (based on the measurements and the scaling 
presented in this study) was found to be even larger than that pre- 
dicted by the transport model. Also the within tree gradients of 
turgor pressure and osmotic concentration can be predicted from 
the axial distribution of xylem sapwood and phloem volumes and 
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FIGURE 8 | Simulated xylem water potential, phloem turgor pressure 
and phloem osmotic concentration axial profiles for cases of phloem 
unloading in sink and phloem unloading along the stem for pine with 
an assumption of maximum sapwood depth of 2 cm (A) and no 
heartwood formation (B). 



area-specific conductivities using a transport model (Figure 8). 
One can hypothesize a feedback loop between local pressure and 
osmotic concentration mediated by xylem and phloem conduc- 
tances, and the local growth rate of new xylem and phloem tissue. 
This feedback loop, spanning several growing seasons, could be 
explained by the direct link between cell division, expansion, and 
cell wall synthesis on the local water and carbon status (e.g., 
De Schepper and Steppe, 2010; Holtta et al., 2010; Pantin et al, 
2012) providing a structural-functional automata for vascular 
development in trees. 

Xylem conductance also decreases with growth in height, but 
not nearly as sharply as phloem conductance. Xylem conduc- 
tance decreases at the rate of square root of tree height, which 
stems from two simple empirical observations: the cross-sectional 
area of xylem sapwood over branching is conserved (Shinozaki 
et al, 1964) and the increase in xylem conduit radius with tree 
height (r x a x 1 ^ 4 ) (West et al, 1999). Various findings also sup- 
port the notion that tree transpiration and photosynthesis rates 
decrease with tree height as either stomatal conductance and/or 
leaf to sapwood area decrease (McDowell et al, 2002a,b, 2005; 
Martinez-Vilalta et al, 2007). The strong decrease of phloem con- 
ducting capacity toward the tree base stems from the observation 
that the phloem width does not change very much from tree top 
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FIGURE 9 | Simulated turgor pressure difference between leaf and root 
phloem as a function of tree height with varying parameterization for 
a case where all phloem unloading occurs in the root (A) and evenly 
along the stem (B) for pine, "sw = 2 cm" refers to a scenario where 
maximum sapwood depth was assumed to be 2 cm, and "no hw" refers to 
a scenario of no heartwood formation, i.e., all of the xylem was assumed to 
be conducting. "K p , tot*2," and " K Pj to t*0.5" refer cases where the absolute 
phloem conductance has been multiplied by a factor of 2 and 0.5, 
respectively. 



to the base (see also Quilho et al., 2000). One possible explana- 
tion for this is that it follows from the limitations of cambium 
activity. Unlike the sapwood, which accumulates over several 
years, phloem apparently needs to be renewed practically yearly 
(Ewers, 1982). Secondary growth results from rate of cell division 
and their subsequent enlargement. These are constrained by the 
length of the growth period and temperature during the period 
but also by the water status and sugar supply to the growth loca- 
tion (e.g., De Schepper and Steppe, 2010; Holtta et al., 2010). 
While it seems that the vigor of the tree may influence the exten- 
sion of growth period (Rathgeber et al., 2011) it is ultimately 
influenced by the duration of the favorable conditions, particu- 
larly in boreal environment. The water status and available sugars 
are influenced by the tree size such that there is less foliage, and 
presumably sugars, relative to stem the bigger the trees. For these 
reasons, the annual width of phloem growth could be limited, 
and although the relatively constant phloem width in axial direc- 
tion means that phloem width increases relative to yearly tree ring 
width from top to base, it is not able to compensate for the dif- 
ferent functional longevities of the tissues. The increasing girth of 
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the trees, however, will help to balance the difference. We could 
even hypothesize that the need to balance the amounts of phloem 
and xylem tissue would be behind the stem diameter growth 
and sapwood turnover. With a given sapwood requirement, its 
higher turnover would necessarily mean faster thickness growth, 
which would have a large impact also on resource allocation 
and tree development (Nikinmaa, 1992). All these functional- 
structural interactions impose strong boundary conditions for 
the tree development and function. 

CONCLUSIONS 

The study showed that important understanding of whole tree 
functions can be gained by dimensional analysis across tree axes. 
The observed regularities in xylem and phloem structure, and 
the Munch- pressure flow hypothesis seem to provide a feasi- 
ble explanation of phloem transport even in the largest trees. 
Sapwood turnover to heartwood seems to have an important 
functional role in affecting the scaling relations for xylem and 
phloem hydraulic conductances and nitrogen allocation. Xylem 
and phloem tissues are clearly a larger sink of nitrogen than the 
foliage as trees grow in height becoming an important and an 
often overlooked factor in the forest nitrogen cycle particularly 
in the nitrogen limited boreal forest where the slow nitrogen 
turnover rate is often the reason for growth limitation. 
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APPENDIX this range. The amount of phloem between 0 and 0.1 m is 

A1. DERIVATION OF THE ANALYTIC SCALING FUNCTIONS FOR XYLEM 

AND PHLOEM PROPERTIES AS A FUNCTION OF TREE SIZE / L \ 5f!l fc m 

The number of furcations at distance of x from the leaf apex is Yf>. 0-Lo — nA p dx =1 — 1 a 2 y 2 L 0 L 0 

(assuming pipe model relations, i.e., nAx is preserved) 



/ i \°pi 

ctzY^ 2 0. 1 m 8fi2 * 0. 1 m (A5) 

\0.1m/ 



= n [e{ A x , ref = n^L^ = L* = 

A x (;t) otiyPijt^ 1 x 5 Pi \x/ Total amount of nitrogen in xylem as a function of tree height is 



Nr. tot = / nA x p N , x p x dx 



where n re f is the number of furcations at the base of tree (n re f = 1) 

is and L is total tree height. *' lul Ji 0 

Total amount of xylem in a tree as a function of tree height L « 

(starting from a distance Lo away from leaf apex) is = / (-) aiY^x^asy^x^pjcdSx 

A 0 W 



Vx, tot 



£nA x dx= 1 a iy ^x^dx = f^L^a^a 3 y^ X -^x^p x dx-. 

J L L^y^x-^dx = f L L^a^dx = f^a^a^x^^dx 



i-L I 8 Pia 1 y f!l c(3y fi3 p x L 8fi3 + 1 

&p 3 + 1 



LWion/'i = cn Y Pl I iPl+1 - ai/»I* 1+1 (A2) = ai t° t3YP 1 3p ^ 8Pl+8P3 + 1 

6th + 1 



aiY Pl «3Y P3 Px 8Pi + 8p 3 + i (A6) 
8fo + 1 0 



If the distance Lois more than zero (e.g., if measurements only 
exist starting from a distance Lrj from the apex as in our mea- 
surements) and if the xylem properties are assumed to remain 

constant between distance 0 and L 0 , the amount xylem between where p x is the density of the xylem tissue (gnr 3 ). We assume p N 

0 and L 0 (0.1 m in this case) must be added to the sum. Xylem to be 0.5 * 10~ 6 gm~ 3 . 

properties were assumed to be constant (the same as at a distance Also, the amount nitrogen in the xylem between 0 and L 0 

Lq from the apex = in this range. The amount of xylem between must be added to the sum. Phloem properties were assumed 

0 and 0.1 m is to be constant in this range. The amount of phloem between 0 

and 1 m is 

V x , o-Ln = nA x dx = ^ a^Lf'Lo / T X Sp 



L °' N Xi 0-L0 = nA xPxPx dx = — oi^Lf'a^pNLf'Lo 

( L \*> V J 

( ) aiy Pl 0.1m 8f!l *0.1m (A3) / T \ spi 

\0.1m) =( ) ' 



Total amount of phloem in tree as a function of tree height is 



Vp, tot 



f L f L /L\ 8 ^ ! 

= / nA p dx = - a 2 Y h x m dx rL f L f L 

JLq JL a \xj N Pi tot = / nA p p p p N dx= / (- 

I I Jo \ x 

-f L^ l a 2 Y h x^ 2 x-^ l dx= f L Spi a 2 y P2 x 8fi2 - 8pi dx = f L 

JL 0 JLa ~ / 

JO 

L Pl «2Y fe f r8fo _ 8gl + 1 r sp 2 -spi + i\ r L XR o o SR , XR XR 

: m _ 8Bl + 1 I 1 - L o J =j L 8 Vy f WW P4+8P2 - 8Pl d* 



aiy pl (l m) 8pl a 3 Y P3 Pivl m* 3 * 1 m (A7) 
Total amount of nitrogen in phloem as a function of tree height is 

L sPl a 2 y te oi4Y p4 PN* 8fl2 *~ 5Pl * 5P4 <fe = 



= _A_M__Mfc-*Pi + l _ i 5p 2 -8p 1 + l\ (A4) _ I s P'a 2 Y P2 «4Pwy P4 Z^ + Sfe-Sp. + l f 8p 4 + sp 2 -»Pi + l\ 

Sp 2 - 8P1 + 1 V 0 / 8p 4 + Sp 2 - &Pi + 1 V 0 / 

(A8) 

where Lo is the height (distance from apex) from which the 

integration is started from. The integration could not be started pjv is the density of the phloem tissue (g m~ 3 ). We assume 

from Lo = 0 as the phloem volume approaches infinity when lo Pjv to be 0.5 * 10~ 6 gm~ 3 , i.e., the same as for the xylem 

approaches 0. tissue. 

Also, the amount of phloem between 0 and Lo must be added Also, the amount nitrogen in the phloem between 0 and Lo 

to the sum. Phloem properties were assumed to be constant in must be added to the sum. Phloem properties were assumed 
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to be constant in this range. The amount of phloem between 0 
and 1 m is 



L x ' ' 1 



N p , o-Lo = TiA p p p p N dx = I — I a 2 Y P2 Lo fc a 4 Y p4 pN4 P4i o 



0.1m 
l*0.1m 



a 2 a80i6Y P2 + 



P8 + 2P(S (0.1m) 8p2+8f!8 + 25p 6 



(A13) 



lm 



a 2 Y P2 (lm) Sfi2 a4Y f!4 pNlm sfi4 *lm (A9) 



Total hydraulic conductance of the xylem is 

K x,\ot = / ( nA xPc,x4) dx 

Jo 

= jf U^j ^ aiY Pl 3c 8Pl a 7 Y p7 x 8te a5Y 2p5 x 28P5 j dx 
= j l (l^a 1 a 7 a 2 5 Y Pl+fi7+2P5 )~ 1 ^ 7 " 28f!5 ^ 



Nitrogen in leaves 

Assuming that leaf area scales with xylem sapwood cross-sectional 
area and that leaf nitrogen content is constant, total amount of 
nitrogen in leaves can be estimated to be 



Aeaf = A x * Q s 



(A14) 



= (L 8 fc ai a 7 a 2 5 yfc + k + 2 fe) ' 



where Q s is the leaf to sapwood ratio 

Meaf. tot = AeafPSLAPN, leaf = -Ax, ref Qs Pleaf PN, leaf (A15) 

In accordance with the pipe model prediction, xylem cross- 
sectional area equals total xylem volume divided by tree height. 
Total leaf nitrogen content can then be written as 



1 /^-jfo-afc+i _ I _ ! p 7 -afc+i\ (A1Q) 

7 - 2Sp 5 + i v 0 / 



-8P 



Also, the xylem resistance between 0 and Lo must be added to the 
sum. Xylem properties were assumed to be constant in this range. 
Xylem conductance between 0 and Lo is 



Meaf, - 



aiV Pi 



QsPSLAPN, leaf (A16) 



-1 



x, 0-i 0 



(«A x p c> x r 2 ) 1 dx 



'a^^+^ + ^i^ + ^ + ^j ifl(A11) 
Total hydraulic conductance of the phloem is 

K p,\ot=f ( nA pPc,pfy dx 

= (l sp i a 2 a 8 a 2 y P2 + p » + 2 P') " ' x~^ + ^ ~ sp » " 2S ^ dx = 



2 + &Pi - 8p 8 - 2Sp 6 + 1 

Also, the amount phloem between 0 and Lo must be added to 
the sum. Phloem properties were assumed to be constant in this 
range. The amount of phloem between 0 and Lq is 



-l 



p, 0-Lo 



("A p p c , P rfj dx 



The parameters used in the calculations are shown in Table 3. The 
values for leaf-to-sapwood area ratio, specific leaf area, and leaf 
nitrogen content were taken for pine from the Supplementary 
material in Martinez- Vilalta et al. (2009) in which these val- 
ues are reported for the SMEAR II Hyytiala research station in 
Hyytiala Southern Finland, the same place from where most of 
the trees were taken for our measurements. Korhonen et al. (2013) 
reported a very similar value, 1.2%, for leaf nitrogen content of 
1.2% at the same site. The corresponding values for aspen could 
not be obtained for the same site, and the values are taken studies 
of Niinemets (1997), and Hoffmann and Usoltsev (2002). Values 
for both pine and aspen are in line with values from metastudies 
containing several species: McDowell et al. (2002a,b) compiled 
data of leaf to sapwood area from 14 different studies contain- 
ing conifer and broadleaved trees. The average of the species was 
~2700 m 2 (leaf area)/m 2 (sapwood). Lusk et al. (2003) compiled 
data average for 23 evergreen angiosperm and 20 conifer popu- 
lations from six sites in American continent and found an SLA 
of 150 g/m 2 and nitrogen content of 1.14% on average. Reich 
et al. (1998) compiled data from 69 species from four functional 
groups (forbs, broad-leafed trees and shrubs, and needle-leafed 
conifers) in six sites in American continent and found an average 
SLA 140 g/m 2 and nitrogen content 1.7%. 

A2. DERIVATION OF THE SCALING FUNCTIONS WITH 
HEARTW00D INCLUDED 

In additional calculations, we also took into account that part of 
the xylem is heartwood. We assumed that there is a maximum 
sapwood radius r slv , max , and that any xylem tissue radial inwards 
to the maximum sapwood depth was heartwood. We assumed 
that heartwood did not contain any nitrogen and did not con- 
tribute to xylem conductance. We did not find analytical solutions 
to the scaling in the case of heartwood inclusion. Results shown 
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are from the following numerical integration, using a numerical 
integration interval (dx) of 1 mm. 



(a lY ^x^ - n( rxyl (x) - r sw , max ) 2 ) a 7 a 5 y^ + 2 ^x h ^ + 2h ^ ' dx 

(A23) 



A sw — A x — Ahy, — otiY^ 1 *^ 1 — Tt (rxyi — r SWt max ) Phloem hydraulic conductance is then 

V 7 (A17) ^ tot = i 0 { nA P ( '- pr p) dx 

The number of furcations at distance of x from the leaf apex is _ f L ( (aiY^ 1 ^ 1 ~ ^(r^\(L) - r SWy max ) 2 ) 

then Jo \ (aiY Pl x 5tSl - it(rxyi(x) - r W| max ) 2 ) 

= "refAsw. ref = («iY Pl ^ - nfrgjgj - r sw , max ) 2 ) c^a^ 2 + p8 + 2p V P2 + m + 28 M " ' <b (A24) 

A sw (x) (aiyPi^Pi - Tt(r xy ,(x) - r sw , max ) 2 ) 7 

Total leaf area is then proportional to sapwood, instead of whole 
Total amount of xylem sapwood is then , , 

' r xylem, cross-sectional area 

v, „ _ f L . I 1 '"■^ i *-r-^- r -"-^> *- - *- * c - ^ 

Jo Jo (aiyPi^Pi - Jt(rwi(x) - r W| max) 2 ) 

A3. DERIVATION OF THE ANALYTIC SCALING FUNCTIONS 

(aiY Pl x 5Pl - nOiyiCx) - r sw max ) 2 ) cfoc FOR THE AXIAL DISTRIBUTION OF XYLEM AND PHLOEM 

PROPERTIES FOR A GIVEN TREE SIZE 



= J (ot lY Pl x sPl - Jt(v(*) - r sw , max ) 2 ) (A19) distribution of xylem volume is 

-1 aiY Pl x 5fil = I 5fSl aiY Pl (A26) 



Vp tot = /" nApdx Note that the distribution of xylem volume is constant with axial 

Jo position. This stems from the pipe model hypothesis. 

,L ( aiY Pi L sPi _ n(r l(L) _ r max) 2) Axial distribution of xylem volume is 

= / 1 (T~mk ; 7S —a 2 Y P2 x° P2 dx 

Jo (aiyf 11 ^ 1 - 1(^1(1) - r W)mal ) 2 ) /L\ 8pi 

(A20) V p (x) = nA p =l-j a 2 y^x m = L^o^ 2 * 8 * 2 " Pl) (A27) 



Total amount of nitrogen in the xylem is then 



Axial distribution of nitrogen in the xylem is 



N x 



X, tot = / tlAxfiN, xPxdX / L\ 1 ;J ;.[.; r. 

Jo V x (x) = nA x p N> x = - n-.^Plv 4 Pln„„P3^P3 



/o VxW = «A x pat, x = I - I oiiY^'x ri c»3Y * 1 



Jo 



( aiy Pi L Spi _ jtCr^iCL) - r 5W , max ) 2 ) = L 5 Pi ai a3Y Pl + P3 * 5p3 (A28) 



0 (aiY pl x 8 Pi - Tt(r xyl (x) - r w> max ) 2 ) , . . , , , 

Axial distribution of nitrogen m the phloem is 

(A21) V p (x) = «A pPN , p = ( - ) 1 oaykx^euY*** 4 



Total amount of nitrogen in the phloem is then = I 5 pi a 2 a 4 Y fc Y P4 * 8 * 2 ~ Pi + ^ ( A2 9) 



f nApppptfdx Axial distribution of xylem conductance is 

Jo 

= f" ( a ^ 1Lm ~ *W> - ^ gg^^yfe^fe p 4 ^ arfx = nA,p c , x r 2 = (-) ' arf 1 &arfi&c& 1 **X a ** 

Jo (a lY Pi^i - n( rxyl (x) - r sw , max ) 2 ) m \x ) 

= I 5Pl a 1 a 7 a 2 Y fSl + P7 + 2fi5 x 5<f!7 + 2P5) (A30) 
Axial distribution of phloem conductance is 

V p (x) = nA pPc , P tj = I - J a 2 Y P2 x 5 P 2 o( 8 Y P8 x 8ps a 2 Y 2fi5 x : 



Np, tot 



(ajyPi^i - n(r^(x) - r sw , max ) 2 ) 

(A22) 

Xylem hydraulic conductance is then 

(nA x p c . x ifj l dx= /rX 5 ^ 



^jc, tot - / ( 
Jo 

= i ((oiiv 1 '^- -tt!/ AVll .vi-,.,, ra , x )-'| = L ; ' p 'a:a,a 2 Y P: + hs + :rv -v' lfl2 - |J, ' +|J,s + ::|Vl (A31) 



L ( (a lV Pl L 8Pl - n(r^(I) - r w , max ) 2 ) 
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